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ON ' 

t-H ' Abstract: 
y—i ■ 

In a previous work, I discussed the effect of the primordial black holes (PBHs) quantum evaporation on 
the reionization of the Universe at small redshifts (z < 60): in principle, the photons emitted during the 
evaporation of such objects could drive a new ionization for the Universe after the recombination epoch 
. (z ~ 1200); this reionization process should happen during the last stages of the PBHs life, when they 

f^S ' totally evaporate and emit a lot of massive and non massive particles. The critical mass of a black hole 

whose lifetime is equal to the present age of the Universe is ~ 4.4 x 10 14 h~ 3 g: thus, PBHs having a mass 
ON ■ M ~ 10 g are the ideal candidates to induce a reionization at small redshifts. 

While in my previous study I considered an exact blackbody photon emission spectrum, here I will 
Oh' adopt a more realistic one, taking into account the quark and gluons jets emission through the contribution 

of a known fragmentation function. When the BH temperature rises above the QCD confinement scale 
Aqcdj one should expect an important contribution from quarks and gluons emission in the form of jets. In 
this paper I also improved my analysis by considering without any approximation the cooling effects in the 
plasma temperature evolution; as a result, I obtained a satisfactory "late and sudden" reionization process, 
, characterised by a very well controlled rise of the plasma temperature: the plasma heating is not so high to 

induce a strong distortion of the CBR spectrum, in agreement with the recent FIRAS upper limit on the 
^V*j , comptonization parameter, y c < 2.5 x 10~ 5 . 
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1. INTRODUCTION 



The possibility that the Universe has been reionized after the recombination is strongly suggested by 
many experimental evidences as, for instance, the Gunn Peterson test [1]: the absence of the Lyman-a 
absorption line in the spectrum of high redshift quasars and the unexpected, low density of neutral hydrogen 
(riHi < 10~ n cm -3 ) in the intergalactic medium (IGM) suggest this hypothesis. 

The causes of such a reionization for the Universe are unclear: the collisional ionization of the IGM by 
cosmic rays [2] or by far-UV photons produced by quasars [3] have been proposed as a possible source but, 
in general, a complete theory describing this phenomenon is not yet available. 

An useful classification of the theoretical models presently known can be found in Ref. [4], where one 
distinguishes "late and sudden" (LS) models (where typically the reionization happens at zr < 60 and it is 
very fast) and "early and gradual" (EG) models, including, for instance, those mechanisms that assume as 
a ionizing source the UV radiation from decaying, massive neutrinos [5]. 

In Ref. [6] I discussed a possible LS reionization mechanism based on the quantum evaporation of 
primordial black holes: as Hawking proved [7], the evolution of a black hole involves a continuous loss of 
mass, a phenomenon ending with the complete evaporation of such an object. This is a typical quantum 
effect: in fact, in a classical sense, a black hole can only absorb particles but the quantum mechanical 
handling of the matter fields in a curved space-time involves some ambiguities in the decomposition of 
the field operators into positive and negative frequency components. Thus, the annihilation and creation 
operators finally work on a vacuum state that is different in regions of the space-time having a different 
curvature [7] and, as a consequence, the creation of particles at the expense of the black hole mass becomes 
possible. 

As I proved in Ref. [6] , the photons emitted during the evaporation process may induce the reionization 
of the Universe; however, in order that such a mechanism works, one should be careful the emitted photons 
do not heat in an excessive way the background plasma. Such a heating should contradict the FIRAS upper 
limit on the comptonization parameter, y c < 2.5 x 10 -5 [8], a value that definitely rules out the possibility of a 
relevant distortion of the CBR blackbody spectrum via the Sunyaev Zel'dovich effect [9]. The proportionality 
of y c to the plasma temperature clearly excludes the possibility that some processes in the past might have 
strongly heated the background plasma: if this is the case, very effective cooling mechanisms should have 
operated to cut off the temperature. 

Due to some numerical problems, in Ref. [6] I solved in an approximate way the coupled differential 
equations system giving the evolution of the ionization degree x and of the plasma temperature T e : thus the 
collisional and excitation cooling terms were underestimated. 

In the following sections, I will present an improved calculation that solves the equation system taking 
into account in an exact way all the cooling terms; moreover, I will consider a modified, non blackbody photon 
emission spectrum which contains also the effect of quantum jets produced by PBHs at high temperatures 
(T > Aqcd)- the importance of such a jets emission has been put in evidence in Ref. [10] and may have 
some influence on the reionization process. 

As a consequence of these improvements, I obtained a satisfactory evolution of both the ionization degree 
x and the plasma temperature T e : in particular, the plasma heating is less relevant than in my previous 
study. The modification of the emission spectrum has the consequence that a lower density of primordial 
black holes is necessary to produce a significant reionization, because in this case the photons coming from 
the secondary decays of hadrons add to the ones produced by the direct evaporation; a rough estimate of 
the present density parameter SIpbh inferred from this calculation is also given. 

Future improvements are still possible: for instance, one may investigate how the production of massive 
particles, in particular electrons, might modify the background plasma and possibly influence the Inverse 
Compton scattering processes. Note however that the emission of charged particles represents a problem 
not immediately solvable because, in this case, one should also study the charge evolution of the black 
hole: from the theory, we know that the loss of charge and angular momentum of a charged, rotating 
black hole is very fast [11] but, as pointed out in Ref. [12], large, electrically charged BHs have a very 
interesting evolution towards the extreme Reissner-Nordstrom limit. During this evolution, very peculiar 
thermodynamic properties appear and various phase transitions may income during the evaporation process. 
The consequences of this particular evolution are a lifetime many order of magnitude larger than the one 
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of a common Schwarzschild black hole and an emission spectrum explicitly containing the charge evolution 
[13]: probably, in these conditions, my analysis of the reionization should be radically modified and the 
conclusions I drew for the uncharged case cannot be simply extended without a deep examination of the 
processes involved in the BHs discharge. 

The consideration of the jets emission represents the first step toward a deeper understanding of the 
problem: in this paper, I will compare the characteristics of the reionization process both in the case of an 
exact Hawking spectrum and in the one when a jet fragmentation function is introduced in order to consider 
also the secondary photon emission from the hadron decays. 

This paper is structured as follows: in Sec. 2, I will discuss the general properties of primordial black 
holes, in particular the characteristic of their quantum evaporation with and without jets emission; in Sec. 
3, I will recall the main formulas and results of my previous work, in particular discussing the differential 
equations system for the variables x and T e ; in Sec. 4 I will present an improved numerical solution method 
that uses the new photon spectrum and does not involve any approximations for the cooling terms; in such 
a way, a noteworthy improvement in the plasma temperature behaviour is obtained. 

An estimate of the density of PBHs one needs in order to have a significant reionization for the Universe 
is also given: it corresponds to a present density parameter Qpbh ranging from 1.12 x 10~ 12 to 1.65 x 1CT 8 , 
depending on the PBHs formation time one assumes. 

Finally, in sec. 5 I will discuss the new results in comparison with the ones obtained in Ref. [6]; the 
conclusion is that the proposed model of a reionization induced by quantum evaporation of PBHs seems to 
work without producing an excessive plasma heating and seems to be effective in the range of reionization 
redshifts [15, 60] here considered. Further extensions to higher zr values are also possible. 

2. THE PHYSICS OF THE PRIMORDIAL BLACK HOLES: 
QUANTUM EVAPORATION AND PHOTON 
EMISSION SPECTRUM. 

The fundamental work of Hawking about black holes physics [7] puts in evidence many interesting 
properties for such objects, in particular their peculiar connections with the general thermodynamics. 

The creation of a black hole becomes possible when a mass contracts to a size less than its gravitational 
radius: typically, this is the case for large mass stars at the end of their evolution; however, for objects 
having smaller masses, a compression to huge densities is necessary in order to create a black hole; in the 
same time, the large pressure forces counteracting the compression should also be overcome. 

For this reason, the formation of black holes having a small mass (M << M ) is very improbable in 
the contemporary Universe: however, Hawking's quantum evaporation phenomena are just important for 
small-mass black holes, the blackbody temperature of the emission being: 



Tic 3 
8irGM 



M 



10 13 g 



-l 



GeV. (2.1) 



Therefore, our interest is mainly fixed upon the so-called primordial black holes (PBHs): at the early stages 
of the cosmological expansion, the density of the matter was really huge, thus enabling the formation of 
small-mass BHs. 

As stressed in Ref. [14], small perturbations in a homogeneous, isotropic, hot Universe would not be 
able to produce relevant inhomogeneities without the contemporary presence of large fluctuations in the 
gravitational field; the possibility of the creation of a black hole just appears when the quantity I = ct 
(t is the time elapsed since the Big Bang) grows to a value in the order of the metric perturbation size 
and its resulting mass will be equal to the mass contained in a volume I 3 at the time t [14]. The role of 
primordial density perturbations in the PBHs formation processes is discussed in a more extensive way in 
[15]; in particular, a study of the effect of the adiabatic inflaton quantum fluctuations in chaotic models can 
be found in [16]; other possible formation mechanisms might be a cosmological phase transition involving 
bubbles collision [17], a softening of the Universe equation of state [18] or a collapse of cosmic strings [19]. 
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In this work, my interest is mainly concentrated on the evaporation of PBHs that form during the first 
stages of the expansion of the Universe: a choice of the possible birth times will be done in the next sections 
in connection with the calculation of the PBHs initial density we need in order to have an appreciable 
reionization. 

The kind of the emitted particles obviously depends on the blackbody temperature of the PBHs: fol- 
lowing eq. (2.1), the mass loss during the quantum evaporation of a BH makes it hotter and hotter, thus 
enabling it to produce more and more massive particles. Typically, BHs having a mass larger than 10 17 g 
emit massless particles only, like photons, neutrinos and, may be, gravitons; electron production should 
be expected by BHs whose mass is in a range 10 15 g < M < 10 17 g and muon production for a range 
10 14 g < M < 10 15 g. 

Here I will consider PBHs that survive till the post-recombination epoch, thus having a mass M <~ 
O(10 14 ) g, near to the critical mass, M c <~ 4.4 x 10 14 h~ 3 g. These PBHs should reionize the Universe at 
a redshift zr that corresponds to the time of their complete evaporation. 

The initial mass Mi is connected to the lifetime of a PBH by the following formula [20] : 



t 



evap 



1.19 x 10 3 



G 2 Mf 
he* /(MO 



6.24 x 1(T 27 fiMi)- 1 Mf sec, 



(2.2) 



where the function f(M) contains the contributions of the different species of particles and it is normalized 
to the unit for very massive (M > 10 17 g) BHs, emitting massless particles only; a suitable way to express 
these contributions is given by the following formula [20] : 
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(2.3) 



1.10- 10 14 

In eq. (2.3) the first addendum in the right-hand side expresses the contribution of electrons, positrons, 
photons and neutrinos; heavier particles contributions are considered in the remaining terms, following their 
relative importance; the factor 3 takes into account the presence of the color charge for quarks and the 
denominators in the exponential terms arc defined as the product (3 S jMj, where Mj is the mass of a black 
hole whose temperature is equal to the rest mass \ij of the j species and f3 S j is a spin-dependent factor 
defined [20] in such a way the energy of a BH having M — f3 S jMj has a peak at fij. 

Eq. (2.3) is fundamental in the determination of the mass evolution of a PBH that evaporates: in fact, 
in Refs. [10], [20], Carr, Mac Gibbon and Webber show that the quarks and gluons emission processes should 
not be neglected at the stages when the BH mass falls below 10 14 g and the temperature T becomes larger 
than the confinement scale Aqcd- at this time one should expect a copious production of jets and this is 
the case we are interested in, due to the particular range of initial masses chosen for the PBHs. 

Through the function f(M), eq. (2.3), we can take into account the quarks/gluons emission processes: 
the mass evolution is given by: 



dM = f 
dt ' ^2nh J 1 
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(2.4) 



eq. (2.4) means that the emission of a parent particle j with total energy Q decreases the BH mass by Q/c 2 ; 
here Tj is the absorption probability for the j particle having a spin Sj [21] and a sum on all the emitted 
species has been performed [20]. After the integration over the energy Q, eq. (2.4) can be rewritten as [20]: 

dM 
dt 



— = -5.34 x 10 25 f(M) M 



~ 2 g sec 1 



(2.5) 
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Now, the Hawking emission rate of particles having an energy in the range (E, E + dE) from a black hole 
having an angular velocity u, an electric potential cf) and a surface gravity n is [7]: 



dN 
~dt 



T dE 
2irh 



exp 



E — nhw — eq 
Hk/2itc 



± 1 



(2.6) 



where the signs ± respectively refer to fermions and bosons and T is the absorption probability of the emitted 
species: for photons, it reads as [21]: 
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(2.7) 



here A is the surface area of the BH and the Planck mass and energy have been introduced in order to work 
with dimensionlcss quantities as in Ref. [21]. 

In most cases, one assumes that the charge and the angular momentum of a black hole are negligible 
because, as Page proved [13], their loss through quantum evaporation happens on a time scale shorter than 
the one characterizing the mass loss; this assumption is certainly reasonable but one should also note that 
in some way it may be relevant to determine the influence of the black hole charge evolution (in particular 
studying the discharge process) in connection with the modification of the electron background; thus, one 
should rather follow the approach suggested in ref. [12]. 

In Ref. [6] I studied the reionization of the Universe induced by quantum evaporation of PBHs by using a 
photon emission spectrum as the one of eq. (2.6), i.e. an exact blackbody spectrum; this assumption enabled 
me just to test if the proposed model of reionization might work at least under very general conditions. After 
this satisfactory test, I will now improve my previous calculation by studying the effect of the presence of 
a jet fragmentation function contribution in the emission spectrum: then, eq. (2.6) should be rewritten as 
[22]: 

dtdE ^ ,/ 2irh \ T / dE ^' K ' 

here x and j respectively refer to the final and the directly emitted particles and the last factor, containing 
the fragmentation function gj x , expresses the number of particles with energy in the range (E, E + dE) 
coming from a jet having an energy equal to Q: namely, it reads as [22]: 



dg JX {Q,E) 
dE 




2m-l 



6(E - km h c 2 ), 



(2.9) 



where rrih is the hadron mass, k is a constant 0(1) and m is an index equal to 1 for mesons and 2 for baryons. 

After determining the dominant contribution to the integral over Q and summing over the final states, 
one can approximate eq. (2.8) by the following formulas [22]: 
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(2.10a) 



(2.106) 



(2.10c) 



where the different expressions refer to the specified Q-value that dominates. A more exhaustive discussion 
of the characteristics of the spectrum given by eqs. (2.10a), (2.10b) and (2.10c) can be found in Ref. [22]. 

In the next Section, I will examine the problem of the reionization of the Universe by considering as a 
source of the ionizing photons the quantum evaporation of PBHs in presence of quarks and gluons jets. 



3. THE TIME EVOLUTION OF THE IONIZATION DEGREE 
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AND PLASMA TEMPERATURE FOR A REIONIZED UNIVERSE. 



As I discussed in Ref. [6], the basic equations that control the ionization degree x and the plasma 
temperature T e evolution with the time are [23]: 

= tpi + t coll ~ t r 1 -> (3-1) 

and 

dT e R T e dx 2 y , 

-jT = - 2 r t < - —)^u + WTx)( r - A ^ ^ 

where t" 1 , t~^ u , t~ x are, respectively, the photoionization, the collisional and the recombination rates, given 
by [23]: 

, 87rm e a 5 , . 1 r max , n 7 

^~-- is --(ir^O(w) 3 <* <") 

tf, _. XU H ,ec- , (1-,, (3.5) 

Eq. (3.1) expresses the combined effects of the photoionization and the recombination and the influence of 
the collisional interactions putting the neutral hydrogen in an excited or ionized state. In eq. (3.2), the 
heating V mainly comes from the photoionization process and reads as 

T = Tpi = 8 ~ 5 m - 1 r max - ( w - A ) 



3V3 



-i rujmax / _ a \ 

(l-x)— dcu [ - 3-i n 7 ; (3.6) 

^* J ujmin ^ 



on the contrary, the cooling A takes into account a lot of contributions having a relevant importance in 
order to limit the rise of the plasma temperature induced by the quantum evaporation of PBHs. These 
contributions are due to the recombination, collisional and excitation processes, Compton scattering and to 
the adiabatic cooling due to the expansion of the Universe. The total cooling reads as 

A = A r + A exc + A co ii + Acompt + A exp , (3.7) 

where, respectively: 

A r = 2 x ID- 1 GeV sec" 1 (^V" (^Y; (3.8) 

\AJ V 0.025/ V 200 / ' 



A. . «r» GeV «- , (. - „ (£A) (££)*, (, 9 > 



A co ;/ = 8 X 10 



Acompt — 10 19 



A„ r = 2. x 1CT 22 GeF sec" 1 (^) (1 + as) (^ )" (3-12) 

(Tcmb is the Cosmic Microwave Background temperature). Both the equations (3.3) and (3.6) explicitly 
contain the photon number density n 7 : the time evolution of this variable is quite difficult to establish 



because it depends on many processes not simply correlated; the best way to describe this evolution is 
through a differential equation that reads as follows [6]: 



drij(u),t) R dn~ ( {uj,t) R n 7 (u),t) 

77 + 75 a ^NOR — 2 — UNOR = 

at R ouj R lo 

/ dn lR dn lPI \ 2 dw dn lPR 

here n 7 pj, n lR are the photon number densities respectively involved in the processes of photoionization 
and recombination and ujnor is a normalization factor, equal to 10~ 6 GeV: this normalization assures one is 
working with dimensionless quantities in the numerical routine that solves the differential equation system. 
The last term in eq. (3.13) is the contribution of the photon source, in our case a number upbh of primordial 
black holes that evaporate. 

From the basic work of Peebles [24], the first term on the right-hand side of eq. (3.14) has the following 
form: 

=C( a e nt - (3 e n ls e-»^ ) - ± (3.14) 
\ at at J at \ J n dt 

Explicitly, a e n 2 e is the rate of recombination to excited states, ignoring the recombination direct to the ground 
state, a e — {crv) is the recombination coefficient and f3 e n\ s e~ ulo '^ T expresses the rate of photoionization; (3 e 
is a constant proportional to a e , n\ s is the population of the Is state of the hydrogen atom, u> a is the energy 
corresponding to the transition L a , T the radiation temperature and n is the nucleon number density. Finally, 
C is a reduction factor which takes into account the effect of the La resonance photons. All the numerical 
values of these quantities can be found in Ref. [25] ; really, we do not mind this contribution because it proves 
to be negligible with respect to the difference n y pi — n lR . 
Now, as in Ref. [6], I put 

n 7 = n lPR - n jPI + n jR , (3.15) 
and with the same substitution performed in [6] I obtain: 

dn 7 (oj,t) n-y r R dui i R 9n 7 (u;,t) 
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a e x n - (3 e n ls e ) - x- + ( -- + co NOR ). (3.16) 

At this point my present analysis of the problem differs from the approach I followed in Ref. [6]: in fact, here 
I will take as a source term in eq. (3.13) the photon emission spectrum corrected by using a jet fragmentation 
function (eqs. (2.8), (2.10a), (2.10b), (2.10c)) while previously I used an exact Hawking photon spectrum, 
eq. (2.6). 



4. THE NUMERICAL SOLUTION OF THE DIFFERENTIAL 
EQUATION SYSTEM. 



I will start the present analysis of the reionization mechanism by evaluating the total photon number 
density coming from the quantum evaporation of PBHs: eqs. (2.10a), (2.10b) and (2.10c) give the emission 
rate for various values of the energy characterizing the emitted particles. 

As I stressed in Ref. [6], I would take into account only the contribution of those photons that may 
really reionize the Universe, i.e., I want to select the components of the emission spectrum really effective by 
comparing the time scale characterizing the photon interactions with the expansion time for the Universe. 
The following condition should be satisfied: 

tint ^ t eX p. (4.1) 

As a consequence of the condition (4.1), the latest time at which the photon interactions may be relevant is 
[10]: 

tfree = ( "X CT ) K (4.2) 
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(nx = number density of the background particles, a = interaction cross-section). 

Among the possible interaction processes, I want to consider the ionization and the recombination: the 
range of energy where these processes are dominating is typically E 1 < 14 KeV but, indeed, the maximum 
allowed energy for photons having a role in the reionization of the Universe is further constricted by the 
condition (4.1); in the pregalactic, matter-dominated era, a photon is affected by ionization losses only before 
a redshift [10]: 

1 + z free ~ 1.0 x 10 19 h- 2 (^) £ 3 5 , (4.3) 

where fl p ~ 0.2 is the present density parameter for protons, h — 0.5 and the energy is expressed in GeV. 

Eq. (4.3) can be inversely used in order to determine E max at a particular value of z, corresponding to 
the reionization redshift zr: in tab. 1 I listed the values I obtained, typically smaller than 8 KeV for a late 
reionization with zr < 60. 

Coming back to the cqs. (2.10a), (2.10b) and (2.10c), due to the upper limit we have for the photon 
energy, the form of the spectrum is the one expressed by eq. (2.10c), holding in the case E << T. 

Taking the jet fragmentation function as in eq. (2.9) and separating the mesons and baryons contribu- 
tions, I can write: 

dn~, , dn~ , dn. 



- - + "Ek\Bar, (4.4) 



diodt ]Tot ~ du;dt lMes ■ dcodt 

(for convenience, in the following discussion the energy of the final photons will be called u). 
For Q = m hadr ~ 300 MeV, eq. (4.4) simply reads: 

<9 n 7 i 1 (, UJ \ 1 / LJ\ 3 

the dominating contribution being the mesonic one. 

Turning now to the solution of the photon density equation, it is useful to adopt dimensionless rescaled 
variables, Co = uj/lj^or with ljnor = 10~ 6 GeV, M — M/10 14 g and t — t/10 15 sec and to write eq. (3.16) 
in the form 

dn 7 2 <9n 7 4 n 7 ^n 1 duj dn^pn I 2 x 10~ 9 dCo \ 
di 3t doj " 3t Cj Cj dt dtdCj y Cj dt J' 

In eq. (4.6) the constant 2 x 10 -9 is a numerical factor coming from the transformation to the new rescaled 
variables and all the negligible terms have been skipped away. 
From the explicit calculation of eq. (4.5), one has: 

dn lPR ^ 2 x 10 15 

dtdCo Cj 1 ■ * 

Now, the general solution of eq. (4.6) can be found in an analytical way by using the method discussed in 
Rcf. [6]: here, I would like just to recall that the general solution of a first order, linear partial differential 
equation that reads as 

Pp + Qq = R, (4.8) 

(P = 7S ' 1 = ff an d ( 3' R are any function of the general variables x, y, z) is a function F(u, v) that 
also solves [26] the equations 

if I = if I = ifl; («) 

here u(x, y, z) = c\, v(x, y, z) = C2 are implicit functions and the constants c\, C2 should be fixed by imposing 
suitable boundary conditions to the particular problem one is studying. 

After the identification of the first variable, x, with the time t, the second one, y, with the energy Cj 
and finally z with the photon number density n 7 , I obtain: 

P = 1 Q = j r (4.10a) 
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the main difference with the expression found in rcf. [6] consists in the modified source photon term. 
Contrarily to my previous analysis, here I do not examine the very particular case when one has equilibrium 
between photoionization and recombination. 

The left-hand side of the equality (4.9) gives: 

^ = (4.11) 
dt 3t 

as in Ref. [6], while for the photon density n 7 I calculate 

3t 



dn-y 3t 1 



4n 7 + — a 



(4.12) 



du 2 u> 
where a — 2 x 10 15 . After solving the linear equation (4.12) one finds: 

3 

n 7 (w) = ——a t + c 2 a) 4 ; (4-13) 
8 

the integration constant c-i is fixed by imposing that the final number of emitted photon is zero for t = t : 

c = I a |, (4.14) 

where to = to/10 15 sec and a good choice for (Do should be done in such a way one has a positive photon 
density at all times; from such a constraint, one obtains iD ~ u! max /12; moreover, a lower limit on the 
photons energy comes from the request that the emission process is not less than the recombination losses: 
this minimum value is equal to 0.20 KeV while, as I discussed in the previous sections, the maximum energy 
depends on the reionization redshift zr. 

I adopted a numerical Merson routine to solve eqs. (3.1), (3.2): note however that, although Mcrson is 
a very powerful method, it fails if the photon number density evolution as expressed by eq. (4.14) is directly 
inserted in in eqs. (3.3), (3.6) for tp] and Tpj . Really, Merson cannot control the contemporary evolution 
of two variables like the time and the photon density in the same equation; a solution is only possible by 
making some simplifying assumptions, namely: 

a) the collision contribution t~ oll in cq. (3.1) is neglected with respect to the photoionization one: this 
assumption is certainly acceptable if we suppose that the quantum evaporation phenomena is the most 
effective contribution to the photon number density. In this way, Merson should not control terms 
having coefficients different by many order of magnitude. 

b) In eqs. (3.3), (3.6) I took Lo min = 0.20 KeV, i.e. the minimum of the energy previously obtained; then 
I put: 

A^t) = P"" dCu ^M; (4.15a) 

Bi{t) =u NO r P"" du ((D - A) (4.156) 

Now, I insert in eqs. (3.1) and (3.2) the time-averaged values A\, B\ that I quoted in tab. 2 for various 
reionization rcdshifts: in such a way, the Merson routine must solve a constant coefficients equation 
system and no approximations are necessary in the handling of the cooling terms of eq. (3.2). 
On the contrary, in ref. [6] I tried to maintain the time dependence in the functions A\, B\ but I was 
obliged to insert some cooling terms in an approximate form; the present choice seems most suitable also 
because the use of an average number of photons in eq. (3.1), (3.2) probably does not represent an extreme 
approximation, due to our poor experimental knowledge of the primeval PBH density parameter Upbh, 
further improvements in the mathematical solution of the system are under study. 
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As a result of this solution method, I am able to control the rise in the plasma temperature T e while 
previously the model involved a too large effect of plasma heating. The effectiveness of the cooling mechanism 
can disclose if the suppression due to the negative exponential terms in eqs. (3.9), (3.10), i.e. to the excitation 
/ collisional terms, is avoided by choosing to start with a high plasma temperature T e corresponding to an 
instantaneous heating. 

Remembering this assumptions, the couple of differential equations is written as follows: 

= C x {z R ) n PBH (1-x) - 77.45 I- T e 1 ; (4.16a) 
at t l 



dT e T e T e dx 0.67 

= -1.33 S- - - + 



C 2 (z R ) npBH (1 - x) - 

/2 



dt ' t (1 + x) di 0- + x) 

- 113.40 t- f e 1/2 - 5.69 x 10 7 x (1 - x) 1 f 1 - ^) 
t 2 t 2 V T e > 

4.54 x 10'. (1 - X) I (1 - ™°) f." 2 - 44.10 X (f. - ^) ^ 

- 2.13 (l + a) -4 



(4.166) 



here T e = Te/lO" 11 GeF. 

In eq. (4.16a) I neglected the collisional term t~ (as I explained above); thus, the addends in cq. 
(4.16a) respectively correspond to , t^ 1 and the ones in brackets in eq. (4.16b) to T pi , A r , A exc , A coU , 
^Compt an d A ea . p . Finally, the coefficients Ci(z R ), C 2 (zr) are given by: 

Ci(z R ) = 7.88 x 10 24 ~Ai{z R ), (4.17a) 

C 2 (z R ) = 7.88 x 10 35 Bi(z R ), (4.176) 

and they are listed in tab. 3 for various z R from 15 to 60. 

In eqs. (4.16a) and (4.16b), the parameter upbh is a weight factor that represents the effective number 
of PBHs in the early Universe we need in order to produce an appreciable reionization effect; a suitable value 
for this free parameter is upbh = 10~ 44 ; an estimate in terms of an effective initial number density of PBHs 
corresponding to this value of upbh can be roughly obtained by considering the following relations: 

< M > npBH 1CA 

" ~ mt-) ■ (4 - 18) 

where < M >^ 10 14 g and 

R{t in ) - Ro (t /t m ) a ; (4.19) 

here Ro = 1.25 x 10 28 cm = 1.4 x 10 10 lyr in a typical cosmological model [27] while for a we can take the 
value —0.5. Then, if one assumes for the PBHs density a behaviour that approximately scales as a power 
with exponent 2/3 of the time and an initial time ti n ~ 10~ 70 -j- 10~ 75 sec, the density parameter flpBH at 
the present time is: 

VLpBH = = 1.12 xlO" 12 ^1.65xl0- 8 ; (4.20) 

Per 

In tab. 4 I listed the results obtained by choosing different formation times for the primordial black holes: this 
rough estimate is purely indicative but in any case it may reproduce quite well the present experimental upper 
limit of the density parameter, derived from the constraints imposed by studying the CMB characteristics 
[10]: n PBH < (7.6 ±2.6) x HT 9 ^-^±0.15). 



5. RESULTS AND CONCLUSIONS. 
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In figs. 1-6 I plotted the behaviour of the ionization degree x vs —z, as results from the solution of eq. 
(4.16a): the six different plots refer to values of the reionization redshift zr ranging from 15 to 60; in figs. 
7-12 I plotted the evolution vs — z of the plasma temperature; for convenience, I also reported in figs. 13-18 
and 19-24 the corresponding results obtained in Ref. [6] by using an exact blackbody spectrum without jets 
emission. 

The present results should be compared with the ones obtained in [6] by a first order solution of the 
differential equation system; we can remark that: 

a) as in Ref. [6], the black holes- induced reionization is only partial; a quite relevant effect is obtained for 
an evaporation redshift corresponding to zr < 30, while for higher values of zr the process of PBHs 
quantum evaporation cannot produce an appreciable phenomenon of reionization; in fact, for zr > 40 
the ionization degree is x < 0.35. However, with respect to the results of Ref. [6], a smaller density of 
PBHs is necessary to reionize the Universe because in this case, one adds to the direct photon emission 
the indirect one, due to the hadrons decays. Moreover, the plasma heating is limited by a powerful 
cooling: the overall effect is that in eq. (3.1) the t~f term is suppressed and the t" 1 term enhanced. 

b) Discussing now the result obtained for the plasma temperature and looking at figs. 7-12, one can remark 
that a real improvement of the calculation has been obtained: the consideration of all the cooling terms 
without any approximations enables me to have a maximum plasma temperature of the order of the 
tenth of eV while previously the PBHs quantum evaporation was producing an excessive plasma heating; 
the irregularities in the plots are due to some numerical fluctuations attributable to the routine solving 
the coupled differential equations system. 

These fluctuations are more evident when a shorter integration time is considered. 

c) Looking at both the results for the ionization degree and the plasma temperature, one can conclude 
that this PBHs-induced reionization mechanism should effectively be "late and sudden" : the possibility 
of an effective reionization for values zr > 30 should be discarded and thus, in this model, it is very 
improbable to have a significative suppression of the CMB temperature fluctuations on small angular 
scales as predicted by Bond and Efstathiou [28] and Vittorio and Silk [29]; in these works, one proves 
that an early reionization could suppress in a significant way these fluctuations, that, following the 
predictions of CDM models and texture scenarios [30], result too large. 

d) The PBHs formation time should be put very far in the past: considering the jet emission contribution 
in the photon spectrum, the best choice of the PBHs formation conditions that enables to obtain a 
well balanced reionization process (i.e. an high ionization degree without an excessive plasma heating) 
suggests a formation time in a range tf orm — 10~ 70 ~ 10~ 75 sec after the Big Bang and an initial density 
Pi = 10 17 -j- 10 24 g/cm 3 (see tab. 4), corresponding to a present density parameter Qpbh ~ 1-12 x 10~ 12 
-7- 1.65 x 10~ 8 . 

Summarizing, in this paper I studied a mechanism of reionization for the Universe induced by the 
quantum evaporation of primordial black holes. Differently from a previous analysis, here an emission 
spectrum taking into account the quark and gluons jets production has been considered. 

A most accurate solution of the system of coupled differential equations, giving the ionization degree and 
the plasma temperature evolution with the time, is performed; as a result, the reionization of the Universe 
happens in an effective, even if partial, way. In the same time, the resulting plasma heating is limited by 
an effective cooling, due to excitation and collisional processes; that avoids to have large (experimentally 
unseen) distortions of the CBR blackbody spectrum via the Sunyaev- Zel'dovich effect. 

As in Ref. [6], I considered here some values of the reionization redshift in the range [15, 60]: the fast 
rise of the ionization degree x seems to suggest that such a model of PBHs-induced reionization should be 
classified as "sudden and late". In fact, the behaviour of x can be approximated by an exponential function: 
this result justifies the analysis of Ref. [31], where I studied the consequences of an exponential reionization 
of the Universe on the polarization of the Cosmic Microwave Background. 

Finally, the possibility to have a total reionization for a redshift zr > 60 seems to be excluded by the 
results here found: in this model, no damping of the temperature fluctuations at small angular scales should 
be expected. 
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Tab. 1: Maximum energy for photons affected by ionization losses before a redshift z. 



z 


Umax (xlO 6 ) GeV 


15 


5.54 


20 


6.00 


30 


6.70 


40 


7.20 


50 


7.70 


60 


8.10 
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Tab. 2: Time averaged values of the integrated photon density A\, B\. 





















ZR 






tR 






Ax 






15 


6.43 


X 


10 15 sec 


10.4 


X 


10 19 


3.83 


x 10 14 


20 


4.30 


X 


10 15 sec 


8.88 


X 


10 19 


3.54 


x 10 14 


30 


2.38 


X 


10 15 sec 


7.12 


X 


10 19 


3.17 


x 10 14 


40 


1.57 


X 


10 15 sec 


6.17 


X 


10 19 


2.96 


x 10 14 


50 


1.13 


X 


10 15 sec 


5.40 


X 


10 19 


2.76 


x 10 14 


60 


0.86 


X 


10 15 sec 


4.88 


X 


10 19 


2.63 


x 10 14 
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Tab. 3: Value of the coefficients C\{zr), C2{zr) respectively appearing in t p \ and Tpj. 



ZR 


Ci 


c 2 


15 


8.20 


xlO 44 


3.02 


xlO 50 


20 


6.99 


xlO 44 


2.79 


xlO 50 


30 


5.61 


xlO 44 


2.50 


xlO 50 


40 


4.86 


xlO 44 


2.33 


xlO 50 


50 


4.25 


xlO 44 


2.18 


xlO 50 


60 


3.84 


xlO 44 


2.07 


xlO 50 



16 



Tab. 4: Values of the density parameter fi pbh and of the PBHs number densities for 
different choices of the initial black holes birth time U (Exp: n PBH < (7.6±2.6) x lCT 9 /^ 1 95±0 15 >). 



U (sec) 


p l (g cm 3 ) 


p a (g cm 3 ) 


QpBH 


10- 75 


4.28 x 10 24 


7.74 x 1(T 38 


1.65 x 1CT 8 


IO- 73 


4.28 x 10 21 


1.66 x 1(T 39 


3.53 x 10" 10 


IO- 71 


4.28 x 10 18 


3.59 x lO" 41 


7.64 x 1(T 12 


IO- 70 


1.35 x 10 17 


5.26 x 1(T 42 


1.12 x 10" 12 
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FIGURE CAPTIONS 



Fig. 1: Evolution of the ionization degree x vs — z induced by primordial black holes evaporation; the Hawking 
emission spectrum is corrected for the presence of quarks and gluons jets; the plot is obtained by 
assuming a reionization rcdshift z = 15. 

Figs. 2-6: The same evolution of x but, respectively, for z — 20, z — 30, z = 40, z — 50, z = 60. 

Fig. 7: The evolution of the plasma temperature T e vs — z; the reionization redshift is z = 15. 

Fig. 8-12: The same evolution of T e but, respectively, for z — 20, z = 30, z — 40, z = 50, z = 60. 

Fig. 13: The evolution of the ionization degree x is calculated by using an exact blackbody emission spectrum in 
a first order recursive calculation (see Ref. 6) and plotted as a function of — z; the reionization rcdshift 
considered is z = 15. 

Figs. 14-18: The same evolution of x but, respectively, for z — 20, z — 30, z — 40, z = 50, z = 60. 

Fig. 19: The evolution of the plasma temperature T e is calculated by using an exact blackbody emission spectrum 
in a first order recursive calculation (see Ref. 6) and plotted as a function of —z; the reionization redshift 
considered is z = 15. 

Figs. 20-24: The same evolution of x but, respectively, for z — 20, z — 30, z — 40, z — 50, z = 60. 
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